schm
schm(
q
b)
Func
Local
q
b1,
q
b2,
q
r,a,ps
ò
(f*g,x,0,1)
à
ps(f,g)
{
q
b[1]}
à
q
b1
{
q
b[1]/(
Ö
(ps(
q
b[1],
q
b[1])))}
à
q
b2
For a,2,dim(
q
b)
q
b[a]-
S
(ps(
q
b[a],
q
b2[j])*
q
b2[j],j,a-1)
à
q
r
augment(
q
b1,{
q
r})
à
q
b1
augment(
q
b2,{
q
r/(
Ö
(ps(
q
r,
q
r)))})
à
q
b2
EndFor
{
q
b1,
q
b2}
EndFunc
